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, A recently introduced lattice model, describing an extended system which exhibits a reentrant 

(symmetry-breaking, second-order) noise-induced nonequilibrium phase transition, is studied under 
the assumption that the multiplicative noise leading to the transition is colored. Within an effective 
Markovian approximation and a mean-field scheme it is found that when the self-correlation time 
r of the noise is different from zero, the transition is also reentrant with respect to the spatial 
coupling D. In other words, at variance with what one expects for equilibrium phase transitions, 
a large enough value of D favors disorder. Moreover, except for a small region in the parameter 
ly-^ ' subspace determined by the noise intensity a and D, an increase in r usually prevents the formation 

ff) . of an ordered state. These effects are supported by numerical simulations. 
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d • I. INTRODUCTION 



In the last few decades we have witnessed a paradigmatic shift regarding the role of fluctuations, from the equilibrium 
picture of merely being a N~ x / 2 perturbation on thermodynamic averages — or triggering at most phase transitions 
between well defined minima of the free energy Q — to lead a host of new and amazing phenomena in far from 
equilibrium situations. As examples we may cite noise-induced unimodal-bimodal transitions in some zero-dimensional 
models 0, stochastic resonance in zero-dimensional and extended systems [&0, noise-induced spatial patterns ||, 
noise-delayed decay of unstable states ||, ratchets 0, shifts in critical points iQ, etc. 

Recently it has been shown that a white and Gaussian multiplicative noise can lead an extended dynamical system 
(fulfilling appropriate conditions) to undergo a phase transition towards an ordered state, characterized by a nonzero 
order parameter and by the breakdown of ergodicity ^ . This result — first obtained within a Curie- Weiss- like mean- 
field approximation, and further extended to consider the simplest correlation function approach — has been confirmed 
through extensive numerical simulations JTo| . In addition to its critical nature as a function of the noise intensity a, the 
newly found noise-induced phase transition has the noteworthy feature of being reentrant: for each value of D above 
a threshold one, the ordered state exists only inside a window [<7i, 02]. At variance with the known case of equilibrium 
order-disorder transitions that are induced (in the simplest lattice models) by the nearest-neighbor coupling constant 
D and rely on the bistability of the local potential, the transition in the case at hand is led by the combined effects of 
D and a through the nonlinearities of the system. Neither the zero-dimensional system (corresponding to the D = 
limit) nor the deterministic one (a = 0) show any transition. 
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This counterintuitive ordering role of noise has also been found afterwards in different models in the literature 
l] M. In Refs. E|,E2[, the authors study a noise- induced reentrant transition in a time-dependent Ginzburg- 



Landau model with both additive and multiplicative noises. Ref. |13| introduces another simple model with a purely 
multiplicative noise, which also presents a noise-induced reentrant transition. This reference also gives evidence that 
the universality class of its critical behavior is that of the multiplicative noise |17 ^] (see also Rcf. Jl9[ for a discussion 
of the universality class of these models). In Ref. fl4j| , a first-order phase transition induced by noise is obtained in 
a system of globally coupled oscillators. A similar first-order phase transition is also found in Ref. |jl5). Finally, Ref. 
|16fl addresses the role of multiplicative noise in the context of phase-separation dynamics. 

Although for the sake of mathematical simplicity all these references (in particular, Refs. PjfOl) studied only the 
white-noise case (the only exception is reference fllqj in which colored noise in space, white in time, is considered) it 
is expected that, because of their nature, fluctuations coupled multiplicatively to the system will show some degree of 
time-correlation or "color" Pj2p|-p3]|, and hence new effects may arise from this fact. For example, a reentrant behavior 
has been found recently as a consequence of color in a noise-induced transition [ p4[ and an ordering nonequilibrium 
pha se transition can be induced in a Ginzburg-Landau model by varying the correlation time of the additive noise 
j25| , p6| . Thus motivated, we have studied the consequences of a finite (but still very short as compared with the 
"deterministic" or coarse-grained time scales) self-correlation time r of the multiplicative noise in systems of this 
kind. We now recall some of the new effects that emerge in this colored-noise case, which have been briefly reported 
in a recent work [E7| : 



Our main finding is that, as a consequence of the multiplicative character of the noise, a strong enough spatial 
coupling D leads invariably (for r > 0) to a disordered state, contrary to what would be expected to happen in 
equilibrium statistical-mechanical models. 

Another important result is that, except for large values of a and very small values of r, color has an inhibiting 
role for ordered states. Moreover, there exists a finite and not very large value of r beyond which order is 
impossible. 

These results represent a major departure from what one can expect on the basis of equilibrium thermodynamics, 
according to which one should tend to think that as D — > oo an ordered situation is favored. Whereas that "intuitive" 
notion remains valid if the multiplicative noise that induces the nonequilibrium ordering phase transition is white 
H[10j, it ceases to be so for r > 0. In the former case, the results could be interpreted in terms of a "freezing" of 
the short-time behavior by a strong enough spatial coupling: for D/a 2 — > oo, the stationary probability distribution 
could be considered to be <5-like, just as the initial one. In our case, an analysis of the short-time behavior of the order 
parameter up to first order in r reveals that the disordering effect of D can only be felt for nonzero self-correlation 
time. As r increases, the minimum value of D required to stabilize the disordered phase decreases and the region 
in parameter space available to the ordered phase shrinks until it vanishes. Thus, the foregoing results can only be 
interpreted once we recall that the ordered phase arises from the cooperation of two factors: the presence of spatial 
coupling and the multiplicative character of the noise (which may eventually lead to "counterintuitive" results). 

It is our purpose in this work to render an explicit account of our calculation and, at the same t ime , to expose and 
to discuss the results more thoroughly. After presenting the model in section ||, we begin section III by introducing 



the approximations needed to render the problem accessible to mathematical analysis. We resort to a mean-field 
approximation like in Refs. |)],[l0| an d to a "unified colored noise approximation" (UCNA) |28|,|29|], devised to deal 



with self-corre lated noises. In section [II A a simplified treatment using the aforementioned approximations is given 



and in section [II B we expose the more sophisticated approach that was actually used to obtain the phase diagrams. 



In this approach the UCNA is complemented with an appropriate interpolation scheme ]29j . In section IV we expose 
and discuss the results obtained within the last approach, comparing the phase diagram with the ones arising from 
the simplified approach and (for r ~ 0) from a perturbative expansion, and the /^-dependence of the order parameter 
m for nonzero r with a numerical simulation |27j . A final discussion of the approach and its results is made in section 



II. THE MODEL 



The model under consideration has been introduced in Refs. |§,|l0 27 1: a d-dimensional extended system of typical 
linear size L is restricted to a hypercubic lattice of N = L d points, whereas time is still regarded as a continuous 
variable. The state of the system at time t is given by the set of stochastic variables {xi(t)} (i = 1, . . . , N) defined at 
the sites of this lattice, which obey a system of coupled ordinary stochastic differential equations (SDE) 
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ii = f(xi) + g(xi)rn + — ^2 ( x i ~ x i) 

jen(i) 



(1) 



(throughout the paper, the Stratonovich interpretation for the SDE's will be meant). Eqs. (|l]) are the discrete version 
of the partial SDE which in the continuum would determine the state of the extended system: we recognize in the 
first two terms the generalization of Langevin's equation for site i to the case of multiplicative noise (j]i is the colored 
multiplicative noise acting on site rj). For the specific example analyzed in Ref. j^], perhaps the simplest one exhibiting 
the transition under analysis (see, however, Ref. Jl3]|), 

f{x) = ~x(l + x 2 ) 2 , g(x) = l + x 2 . (2) 

The last term in Eqs.(^) is nothing but the lattice version of the Laplacian \7 2 x of the extended stochastic variable 
x{v, t) in a reaction-diffusion scheme. n(i) stands for the set of 2d sites which form the immediate neighborhood of 
the site r,, and the coupling constant D between neighboring lattice sites is the diffusion coefficient. 

As previously stated, it is our aim in this work to investigate the effects of the self-correlation time r of the 
multiplicative noise on the model system just described. To that end we must assume a specific form for the noises {r]i}: 
we choose Ornstein-Uhlenbeck noises, i.e. Gaussian distributed stochastic variables with zero mean and exponentially 
decaying correlations: 

(Vi(t)Vj(t')) = Ma 2 /2T)exp(-|t-;'|/T)- (3) 

They arise as solutions of an uncoupled set of Langevin SDE: 

rrji = -r\i + cr& (4) 

where the {£,i(t)} are white noises — namely, Gaussian stochastic variables with zero mean and 5-correlated: 
(£i = SijS(t — t'), For r — > 0, the Ornstein-Uhlenbeck noise T)i(t) approaches the white-noise limit (t) with 
correlations (t)^ (t')) = a 2 5ij S(t — t'), which was the case studied in Refs. 



III. THE APPROXIMATIONS 



The non-Markovian character of the process {xi} (due to the presence of the colored noise {rji}) makes it difficult 
to study. However, some approximations have been devised which render a Markovian process that — whereas nicely 
simplifying the treatment — can still capture some of the essential features of the complete non-Markovian one. The 
aforementioned UCNA is one of them: in fact a very reliable one, because of its ability to reproduce the small- and 
large-r limit s |28| . By resorting to interpolation schemes, one can retrieve meaningful results in wider vicinities of 
these limits p9|. 

As already declared, our approach is a mean-field-like one. The earlier we make this kind of assumptions in the 
calculation, the cruder the approximation will turn out to be. In order to find the phase diagram in the presence 
of colored noise we have made the mean-field approximation at some late stage, so enhancing the precision of the 
calculation. However, since this calculation is a tedious one, we shall first expose a simpler approximation which 
brings out most qualitative results. We aim in this way to underline the physical origin of the results presented in 
section IV. The differences arising from both calculations are pointed out there. 



A. A simpler approach 

The simpler mean-field approximation follows closely Curie- Weiss' mean-field approach to magnetism, and consists 
in replacing the last term in Eqs.(Q) 

Al ~ 2d ^ ( Xj ~ ^ 

jen(i) 

by 

Ai = D(x-Xi), (6) 
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where x is a parameter that will be determined self-consistently. In other words, the (short-ranged) interactions are 
substituted by a time- and space-independent "external" field whose value depends on the state of the system. Since 
in this approximation Eqs.([iJ) get immediately decoupled, there is no use in keeping the subindex i and we may refer 
to the systems in Eqs. ([j]) and (Q) as if they were single equations. Hereafter, the primes will indicate derivatives with 
respect to x (clearly A' = —D). 

If we take the time derivative of Eq.(Q), replace first i) in terms of n and £ from Eq.(Q) and then r\ in terms of x 
and x from Eq.(Q), we obtain the following non-Markovian SDE: 



T(X - ^X 2 ) = — [l — T 

9 



(f + A)' - y -(.f + A) 
9 



+ (f + A) + o-gii. 



(7) 



The aforementioned UCNA allows us to recover a Markovian SDE: for our particular problem it amounts, on one 
hand, to a usual adiabatic elimination (namely, neglecting x) and, on the other, to neglect x 2 so that the system's 
dynamics be governed by a Fokker-Planck equation [^0[ . The resulting equation, being linear in x (but of course not 
in x), can be immediately solved for x giving 



with 



x = Q(x; x) + S(x; x)£, 



Q(x;x) = (f + A)9, 
S(x; x) = agO, 

e(x;x) = {l-rg[(f + A)/g]'}- 1 - 



The parametric dependence of Q(x) and S(x) on x has been written explicitly. 
The Fokker-Planck equation associated to the SDE Eq. (ph is 



d t P(x,t;x) 

with drift and diffusion coefficients given by 



d x [R\(x\ x)P(x, t; x)} + -d 2 [i? 2 (^; x)P(x, t; x)] , 



f 



R 1 (x;x)=Q + -(S 2 y 
R%(x; x) = S 2 . 

The solution of the time-independent Fokker-Planck equation leads to the stationary probability density 

^R^x'-x) -d x .R 2 (x';xY 



P st (x;x)=JV- 1 exp 



dx' 



R2(x'; x) 



(8) 



(9) 
(10) 
(11) 



(12) 

(13) 
(14) 

(15) 



being AT its norm. The partial-derivative notation d x i in Eq.(p"5|) — as in Eqs. ([l8]) and (|T^) below — is only a reminder 
of the parametric dependence of R\, R2 on x. 

The value of x arises from a self-consistency relation, once we equate it to the average value of the random variable 
Xi in the stationary state 



(x) 



dxxP st (x;x) = F(x). 



As in the known Curie- Weiss mean-field approach, the condition 

= 1 



dF 
dx 



(16) 



(17) 



allows us to find the transition line between the ordered and the disordered phases. Here also, F(x) is a smooth odd 
function such that Eq.([l6|) has always a root at x = and for dF/dx\ s _ > 1 it has two nontrivial roots which differ 
only in sign. The condition Eq.([j"7j) thus reads: 



TV" 



dx x dx' exp 



dx 



11 2i?i — d x nR 2 



Ri 



Of 



{ 2Ri — d x <R 2 



R-2 



(18) 
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where 



(19) 

2=0 

Eqs. (PI) and (|l^) must be solved numerically in order to find the lines in parameter space (<t,t,D) that separate 
ordered (x ^ 0) from disordered (x = 0) phases. The results of this calculation will be shown in section V. Next, we 
introduce a more refined approach in which the mean-field approximation is made at a later stage in the calculation. 



Af ■■ 



dx exp 



dx 



, 2Ri — d x 'R2 



B. A more refined approach 



As we shall see, the relations obtained with this more sophisticated approach are similar to Eqs. (|T^) through JlS|), 
but with different expressions for the functions Ri(x; x) and R2(x; x). The idea here is to introduce first the UCNA, 
without yet resorting to the mean- field approximation. In the following, A, has the same meaning as in Eq.(^) and, 
as it occurred previously with A', it satisfies A^ = —D. Note however that whereas for /, = f(xi) and gi = g(xi) 
the prime has the meaning of a total derivative with respect to Xi, for A^ and all the functions involving them, its 
meaning is really that of a partial derivative with respect to a;,. Proceeding as before, i.e. taking the time derivative 
of Eqs.ffl) and using Eqs. (Q) and to eliminate the ry's in favor of the x's and £'s, we obtain the following system 
of (non Markovian) SDE's 



t(x< - ^xf) = 
9i 



1 - TQi 



ft + A, 
9i 



ii + (fi + Aj) + agi£i + — ^2 

j'6n(») 



(20) 



The UCNA proceeds here through the neglect of xi and of (ii) 2 , so retrieving a linear equation in the x's (but of 
course not in the x's), which can be rewritten as 



[Qi + Sib] 



fi 



{fi + A,) + ag^i + 



Dt 
~2d 



*: 

j'en(i) 



ft 

2d 



j£n(i) 



Here the quantities 



9i = [i - T 9i (fi + Ai/ gi y 

Qi = (fi + Ai)0i, 
Si = crgiOi, 
ft = DrOi 



(21) 



(22) 
(23) 
(24) 
(25) 



have been defined in order to simplify the notation. Note that although only the dependence upon Xi has been made 
explicit in this notation, these quantities depend also (through A^) on the values Xj at the neighboring sites. 

Now, assuming the lattice to be isotropic, we apply to this set a mean-field-like approximation (but not yet the 
main one) consisting in replacing in all the functions appearing in Eq. ( ^l|) the 2d neighbors Xj of the variable x% by a 
common value y^. Through this procedure one reduces to two the number of different coupled SDE's: one for x = X4 
and another for its nearest neighbor variable y = j/,-. These are 



a = h a + g ab £ fc 



(26) 



where a sum over the values x, y is implied for the indices a, b, and the noise variables satisfy (£, a (t)£,b(t')) = <>ab 5(t — t'). 
If, similarly as before, we define 



9(x,y) 



d f f(x)+D(y-x) 
'-^dx-y gjx) 



Q(x,y) = [f(x)+D(y-x)]6(x,y), 
S(x,y) = ag(x)9(x,y), 
P(x,y)=TD9(x,y), 

A(x, y) = [1 - (3(x, y)f3(y, x)}' 1 = A(y, x), 



(27) 

(28) 
(29) 
(30) 
(31) 



5 



and write a = y if a = x and viceversa, then the explicit forms of the functions in Eq. (|26|) are 

h a = A(x, y) [Q(a, a) + 0(a, a)Q{a, a)] (32) 

and 

g ab = A(x,y)S(a,a) if b = a, (33) 

= A{x, y) /3(a, a)5(a, a) if b = a. (34) 

The bivariate Fokker-Planck equation associated to Eqs.(p6|) is 



d t P = -d a (R a P) + ~d a d b (R ab P), (35) 

where P = P(x, y; t). According to standard techniques |0J, the drift and diffusion coefficients are given respectively 
by 

R a (x,y) = h a + ^gbcdb(gac), (36) 
R ab (x,y) = g ac gbc- (37) 
Since the denominators occurring in these equations may become zero for some values of x or y, we resort to an 



interpolation procedure (analogous to the one used in Refs. [|24 29 1) consisting in replacing the expression (J3l|) for 
A(x, y) by 

_ l-/3{x,y)0(y,x) 
A{X > y) -l + P(x,yrt3(y,xr (38) 



Since f3{x, y) is proportional to r, it follows that the expression in Eq.(|3l|) coincides with the interpolated one, Eq. 
as t — » and r — > oo (the latter limit meaning indeed "r comparable with the 'deterministic' time scales"). 

By integrating the bivariate Fokker-Planck Eq.(|35|) with respect to y we obtain a single- variable equation, which 
in the stationary case reads 

= -d x (R 1 P st ) + ±d 2 x (R 2 P st ), (39) 
being RiP st and R 2 P st functions of x only: 

/oo 
dyP st (x,y)R x (x,y), (40) 
- OO 

/CO 
dyP st (x,y)R xx (x,y). (41) 
-oo 

Here it is when we resort to the main mean-field-type approximation, resembling the Curie- Weiss' type of approach 
used in Ref. M: assuming P st (x,y) w P st (x)8{y — x), R\ and R2 in Eqs. (|4(]) and ( pil| ) are approximated by 

Rx =R x (x;x), (42) 
R 2 = R xx (x;x). (43) 

In this way, from the stationary joint probability density function P st (x,y) we retrieve an effective single- variable 
one P st (x; x) whose expression in terms of R\(x; x) and R 2 {x\ x) arising from Eq.(|39|) is the same as in Eq. (|l5|) . The 
value of x follows again from a self-consistency relation like Eq. ([l6]) . The procedure to find the phase diagram is the 
same as in the foregoing subsection and the explicit expression of the condition dF[dx\ s=Q = 1 is given by Eqs. jl^ ) 
and (|l9|), this time in terms of the corresponding Ri, R 2 given by Eqs. (42) and JH). 

As discussed in Ref. |27|, although the kind of approximation leading to Eq.(|3§|) is not of a perturbative nature, 
it has provided sound results in the cases analyzed heretofore [ p4|j29| ]. Nonetheless, for the sake of comparison we 
have also adapted to the case of a multiplicative noise a known perturbative procedure p3|. Within this context, the 
expressions for R\(x\x) and R 2 (x\x) come from a consistent small-r expansion of a Fokker-Planck equation (Eqs. 
©or®): 

Ri(x; x) = (f + A) + a 2 x{g[l + r(f + A)'] + r(f + A)}, (44) 
R 2 (x;x) = {<jg[l + T(f + A)'}} 2 . (45) 



G 



IV. THE RESULTS 



A. Phase diagram 

In the following we shall describe the results obtained through the numerical solution of Eqs. ( fL8| ) and (|l^) in the 
more refined approach, i.e. with R\ and R2 as prescribed by Eqs. (ji^) and (^3|). We shall also compare these results 
with the ones arising from Eqs. (13) and (14), and with a perturbative expansion for small r. Figures [l] to || are 



respectively the projections onto the a-D, r-a and t-D planes, of cuts of the boundary separating the ordered and 
disordered phases performed at fixed values of the remaining parameters. 

Let us first concentrate on Fig.[l]: it corresponds to Fig.l in Ref . Jgjjl , but is the result of an improved calculation 
based on the more refined mean-field approach described in Sect. pHB . The novelty is that, at least for r not too 



small, it is now evident that the region available to the ordered phase is bounded. The noteworthy aspects are the 
following: 

A. As in the white- noise case r = (Rcfs. P,|lO||), the ordering phase transition is reentrant with respect to a: for 
a range of values of D that depends on r, ordered states can only exist within a window [a\, a^\. The fact that 
this window shifts to the right for small r means that, for fixed D, color destroys order just above o~\ but creates 
it just above a 2 - 

B. For fixed a > 1 and r ^ 0, ordered states exist only within a window of values for D. Thus the ordering phase 
transition is also reentrant with respect to D. For r small enough the maximum value of D compatible with the 
ordered phase increases rather steeply with a, reaching a maximum around a ~ 5 and then decreases gently. 
For t > 0.1 it becomes evident (in the ranges of D and a analyzed) that the region sustaining the ordered phase 
is closed, and shrinks to a point for a value slightly larger than r = 0.123. 

C. For fixed values of a > 1 and D larger than its minimum for r = 0, the system always becomes disordered for r 
large enough. The maximum value of r consistent with order altogether corresponds to a ~ 5 and D ~ 32. In 
other words, ordering is possible only if the multiplicative noise inducing it has short memory. 

D. The fact that the region sustaining the ordered phase finally shrinks to a point means that even for that small 
region in the a-D plane for which order is induced by color, a further increase in r destroys it. In other words, 
the phase transition is also reentrant with respect to r. For D large enough there may exist even two such 
windows. 

Some of the features just described become more evident by looking at Fig.|[ 

• the existence of a maximum correlation time consistent with ordering for each value of D (occurring for an 
optimal value of a) (Q ; 

• the ordering ability (^) of a very small amount of color for a > o~2(D) (02(1", D) increases very rapidly at first); 

• the reentrance with respect to r and even the occurrence of a double reentrance for D large enough Q). 

Fig.|| represents another way of seeing the reentrance with respect to D for constant a (large enough) and the fact 
that there exists a maximum r consistent with order for each value of a (being it largest for a ~ 5). The scarce 
dependence of D on r in the lower branch — as well as its almost linear dependence on a — is easily understood by 
looking at the rightmost branch of Fig.[l]. 

In Figs . ^ to ^ we compare the results just shown — obtain ed as we said using the the more refined approach of 
Sect JlII B — with the ones arising from the simpler one (Sect. [nT"A ). Figure || corresponds to Fig.§, whereas Figs. 



H and |6| focus respectively on the r = 0.03 and t = 0.015 curves in Fig.|l|. Not only does the simpler approach 
(grossly) overestimate the size of the ordered region but also — as one may infer from Figs. f| and g — it seems to 
predict unbounded ordered regions. 

Figure p| co rresponds to a rather small value of r, so that a comparison with the results obtained by using expressions 
( p4| ) and p5| ) makes sense. For a and D small enough the curves almost coincide. As it is well known, the simultaneous 
consideration of small t and large a cannot be done independently [2l|-|23]|. In the present case, a similar effect arises 
when we consider large values of D, as we discuss below. Hence, in order to consider larger values of a and D, one 
should take extremely low values of r. As Fig.^ shows, already for r = 0.015 there is an apparent discrepancy between 
the perturbative results and the mean-field ones even for not so large values of a and D. The noteworthy fact is that 
the perturbative expansion also tends to indicate the existence of a reentrance with respect to D. 
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B. Order parameter 



The order parameter in this system is to = |x| namely, the positive solution of the consistency equation (see Eq.(|lq) 



in Sect. [II A ). In Figj7j we plot to vs. a for D = 20 and different values of r. Consistently with what has been 
discussed in (Q) and (Q and shown in Fig.|l], we see that as t increases the window of a values where ordering occurs 
shrinks until it disappears. One also notices that at least for this D, the value of a corrresponding to the maximum 
order parameter varies very little with r. Figure ^ is a plot of to vs. r for D = 45 and two values of a (~ 7.07 and 
w 8.94) that illustrates the case of double reentrance in r. 

Since the previous results have been obtained in the mean-field approximation, we have also performed numerical 
simulations in order to have an independent check of the predictions. As a representative example — corresponding 
to the phenomenon (|^) above (the destruction of the ordered phase by an increasing coupling constant D) — we plot 
jointly in Fig|| for rather small values of a (= 2) and r (= 0.01), the D-dependence of the order parameter as 
predicted by our mean-field approach and as resulting from a numerical integration of the original SDE's, Eqs. ([!]). 
We have taken three different lattice sizes in order to assess finite-size effects. As we see, the numerical simulations 
do predict the disordering for large enough D, and even the maximum ordering occurs in a region which is consistent 
with the mean-field prediction. This comparison warns us, however, that the mean-field approximation can severely 
underestimate the size of the ordered region. 



The short-time evolution of (x) can be obtained multiplying Eq.(12) by x and integrating: 

d(x) 



dt 



dxR 1 {x;x)P{x,t;x). (46) 



Let us assume an initial condition such that at early times P(x, t ~ 0; x) — 5(x — x). Equating x — (x) as before, we 
obtain 

^ = Ri(x,x) (47) 

(again, we can use for R\ the result Eq.([l3|) of the simple approximation or the more elaborate one given by Eq. (fL?|) ) . 
The numerical solution of Eq.(|4^) has an initial rising period (it is initially unstable) reaching very soon a maximum 
and tending to zero afterwards. 

For D/a 2 — > oo, Eq. (|47j) results to be valid also in the asymptotic regime since P st {x) = 8(x — x) |10|. In Ref. E^] 



an equivalent equation is obtained in the r = case for both limits (D = and D/a 2 — > oo) being there interpreted 
in terms of a "freezing" of the short-time behavior. According to this criterion, in the D/a 2 — > oo limit the system 
undergoes a second-order phase transition if the corresponding zero-dimensional model presents a linear instability 
in its short-time dynamics, i.e. if after linearizing Eq.(p^): 

(x) = -a(x) (48) 

one finds that a < 0. We then see that the trivial (disordered) solution (x) = is stable only for a > 0. For a < 
other stable solutions with (x) ^ appear, and the system develops order through a genuine phase transition. In 
this case, (x) can be regarded as the order parameter. In the white noise limit r = this is known to be the case for 
sufficiently large values of the coupling D and for a window of values for the noise amplitude a <E [ax, 02]. 

We discuss now how the stability of the ordered phase is altered by nonzero values of r. If we linearize Eq.(|4^) 
using the expression of i?i(x,x) from Eq.(O), we obtain 



{1 + t + tD) 3 

If we use instead Eq.(f42|), the result can be written exclusively in terms of tD 



(l + r + TD) 2 -a 2 (l-3r + 2rD) 



*B(tD) 



with 



A(x) = 1 + 5x + 8x 2 + 3x 3 - 3x 4 - x 5 + 2x 6 + x 7 , (51) 
B(x) = 1 + 3x + x 2 - 5x 3 - 3x 4 + 3x 5 + x 6 , (52) 
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so the instability occurs at a 2 = A/B. Now, the ratio B/A keeps always below 1 in the positive range, has a minimum 
of ~ 0.05 at tD ~ 1.09 and a maximum of ~ 0.36 at tD ~ 2.33. If we consider e.g. the conditions in Fig.|| (namely 
t = 0.03) and take D large enough so that tD > 2.33, one can see that Eq.(|50|) approximates the left boundary better 



than Eq.(49) does. In the limit r <C 1 (but still finite) Eq.(|50|) can be approximated to the expression reported in 



Rcf. [g7|, namely 

1 + tD - a 2 



1 + tD 



(53) 



It is worthwile to stress the fact (evident from Eqs. (^), (^0|) and ([33|)) that the value of D has no effect on the 
location of the instability by itself, but only through the combination tD: according to Eq.(p3|) the stable phase is the 
disordered one ((x) = 0) for 1 + tD > a 2 (since a > 0) and the ordered one ((x) ^ 0) for 1 + tD < a 2 . In summary, 
whereas the noise intensity a has a destabilizing effect on the disordered phase, as soon as r ^ the spatial coupling 
D tends to stabilize it. For r = the last effect is not present, being then a > 1 the condition for ordering P,[To||. 
Considering that the effect of even a tiny correlation is enhanced by D, we can understand the abrupt change in slope 
(from negative to positive) shown in Fig.[l] as soon as r ^ 0. Note the approximately inverse relation between r and 
D for fixed a on the upper branches of Fig.[| even when Eq.(|I§|) is strictly valid for D/a 2 — » oo. 



V. CONCLUSIONS 



This work has focused on the effects of a self-correlation in the multiplicative noise, on the reentrant noise-induced 
phase transition introduced in Rcf. |)| . Whereas in a recent Letter we have reported the most relevant results [^7| , it 
has been our purpose in the present work to expose in more detail the techniques and the approximations employed. 
We also discuss more thoroughly the results, adding new figures and enriching the contents of others. 

Through the use of the UCNA we recovered a Markovian behavior for the system, and through an interpolation 
scheme similar to the one introduced in Ref. ]29| ] we resolved indeterminacies in the equations describing it. We stress 
that the equation resulting from this interpolation scheme are exact in the limits t = and r — > oo. In addition 
to the fact that the interpolation scheme has been already applied with success in other works |24|,|2{|, the goodness 
of our approximations for small but nonetheless finite values of r has been checked against a standard perturbative 
expansion ]32| (adapted for multiplicative noise). It is worth to emphasize the fact that these approximations are, so 
far, the only tool available for an analytical treatment of this essentially non-Markovian problem. 

The main result is that for t / 0, the order established as a consequence of the multiplicative character of the noise 
can be destroyed by a strong enough spatial coupling. Figure [l] shows that for given r (0.03) and <r, the ordered phase 
can only exist between definite values of D. In particular, the upper bound on D decreases roughly as r _1 for given 
a. 

The foregoing result can be understood by recalling that the ordered phase arises as a consequence of the collabo- 
ration between the multiplicative character of the noise and the presence of spatial coupling. The disordering effect 
of D arises only when r ^ (the results in Ref. — rightly interpreted in terms of a "freezing" of the short-time 
behavior by a strong enough spatial coupling — are thus consistent with ours). As r increases, the minimum value 
of D required to stabilize the disordered phase decreases rapidly, and the region in parameter space available to the 
ordered phase shrinks until it disappears. 

The example worked throughout this paper shows vividly the fact that the conceptual inheritance from equilibrium 
thermodynamics (though often useful) is not always applicable. The equilibrium-thcrmodynamic lore would induce us 
to think that as D — > oo an ordered situation is favored |l| . Although this is certainly true for the Curie- Weiss model 
(since in that case the deterministic potential is itself bistable and an increase of spatial coupling has the effect of 
rising the potential barrier between the stable states), it is not in the case we are dealing with, since the deterministic 
potential is monostable. Hence, it is the combined effects of the multiplicative noise and the spatial coupling that 
induce the transition. 

As a summary, whereas one might say that the value of Refs. |^,|T2|,|l0j is that they tell experimentalists where 
not to look for a noise-induced phase transition — namely, in those systems which are prone to exhibit a usual (zero- 
dimensional) noise-induced transition, and for too large noise intensity — the present work tells moreover that, due to 
the consideration of the more realistic colored noise source, an ordered phase is not to be found for large values of the 
spatial coupling either. Though the specific choice of the forms for the functions f(x) and g(x) may appear to some as 
physically unmotivated, is up to our knowledge the simplest one exhibiting this phenomenon. We believe nonetheless 
that the phenomenon is robust and transcends the specific choice made in this work. 

The next obvious step is to consider a finite correlation length in the lattice model, which requires to go beyond 
the mean field approach. This problem is being presently studied. 
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FIG. 1. Phase diagram in the a-D plane, for different values of r: (1) r = 0; (2) r = 0.015; (3) r = 0.03; (4) r = 0.05; (5) 
r = 0.1; (6) t — 0.123. The ordered phase exists only inside the corresponding curves. 
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FIG. 2. Phase diagram in the r-a plane, for different values of D: (1) D = 20; (2) D = 32; (3) D = 45. The ordered phase 
exists only inside the corresponding curves. 




FIG. 3. Phase diagram in the t-D plane, for different values of a: (1) a 2 = 10; (2) a 2 = 20; (3) a 2 = 30; (4) a 2 = 50. The 
ordered phase exists only inside the corresponding curves. 




FIG. 4. Comparison between the simpler mean-field approach (solid line) and the refined one (dots) in the r-cr plane, for 
D = 20 (lower two curves) and D — 45 (upper two ones) . Not only does the simpler approach tend to overestimate the size of 
the ordered region, but it even predicts unbounded ordered regions for some values of D. 
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FIG. 5. Comparison between the two mean-field approaches in the a-D plane, for r = 0.03. Solid line: simpler one; dots: 
refined one. 
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FIG. 6. Comparison between the two mean-field approaches and with a perturbative expansion in the a-D plane, for 
r = 0.015. Solid line: simpler MF; dotted line: refined MF; dashed line: perturbative. 
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FIG. 7. Order parameter m vs. cr, for D — 20 and four values of r: (1) r = 0.015; (2) r = 0.05; (3) r = 0.06; (4) r = 0.07. 
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FIG. 8. Plot of m vs. r for D — 45, showing cases of double reentrance: (1) a 2 = 50; (2) a 2 — 80. 
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FIG. 9. Plot of m vs. D for <j = 2 and r = 0.01, showing the predictions of the more refined mean-field approach together 
with results coming from a numerical integration of the original SDE's, Eqs.(fil). 



15 



